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The responses of the electrical current and the current-induced spin transfer torque (CISTT) 
to an ac bias in addition to a dc bias in a magnetic tunnel junction are investigated by means 
of the time-dependent nonquilibrium Green function technique. The time-averaged current (time- 
C ' averaged CISTT) is formulated in the form of a summation of dc current (dc CISTT) multiplied by 

3 . products of Bessel functions with the energy levels shifted by mhcoo. The tunneling current can be 

viewed as to happen between the photonic sidebands of the two ferromagnets. The electrons can 
pass through the barrier easily under high frequencies but difficultly under low frequencies. The 
tunnel magnetoresistance almost does not vary with an ac field. It is found that the spin transfer 
torque, still being proportional to the electrical current under an ac bias, can be changed by varying 
frequency. Low frequencies could yield a rapid decrease of the spin transfer torque, while a large 
ac signal leads to both decrease of the electrical current and the spin torque. If only an ac bias is 
present, the spin transfer torque is sharply enhanced at the particular amplitude and frequency of 
the ac bias. A nearly linear relation between such an amplitude and frequency is observed. 



PACS numbers: PACS numbers: 73.40.Gk, 73.40.Rw, 75.70.Cn 



C3 . I. INTRODUCTION 



The spin-polarized electrical transport in magnetic multilayer structures has received much attention both exper- 
imentally and theoretically in the last several years [J. The yield of the extensive investigation promises fascinating 
Q implication for applications in information technology. For instance, people may make use of the different resistive 
O states corresponding to parallel and antiparallel magnetizations of different layers as memory elements in magnetic 
random-access memories (MRAM) . It is now known that the antiparallel alignments of the moments in magnetic layers 
I ■ lead to a higher electrical resistance than the parallel alignments, giving rise to a so-called giant magnetoresistance 
J>. , (GMR) or tunnel magnetoresistance (TMR) depending on different magnetic junction systems. The origin of these 
f"~^ ■ phenomena is widely believed to be mainly caused by the spin-dependent scatterings of the conduction electrons. As 
l/~) ' is pointed out '^, when the junction bias is increased, the phenomenon is significantly reduced. Generally speaking, 
Cn the resistive state can be changed by applying an external magnetic field, because the latter could lead to variations 

of the traversing paths of the conduction electrons. However, it is now accepted that not only the electrical current 
is strongly affected by a magnetic state, but also the electrical current can conversely control the magnetic state of 
^—s ■ the magnetic junctions. This effect is predicted by Slonczewski ^ and Berger |j| in magnetic multilayer systems by 
■~»^ , noting that a spin-polarized electrical current can transfer local spin angular momenta of incident electrons to the 
C^ • scattering ferromagnet, thereby exerting a torque on the magnetic moments and therefore changing the magnetic 
^ [ state. When the current is large enough, the magnetization of ferromagnetic layer can be switched. As a result, the 
I spin transfer effect may provide a mechanism for a current-controlled magnetic memory element. Experiments have 

' O ^ given evidence of this effect in the Cu/Co multilayers yl, nickel nanowires la|, manganite junctions i2l]j point contact 
C , magnetic multilayer devicesjg, bulk manganese oxidesQ, tunnel junctionsQ], and Co/Cu/Co spin valve devices [lO| . 
etc. A thermally activated switching of magnetic domains in a Co/Cu/Co spin valve is also confirmed 'lll|. 

To deal with the spin transfer effect, it is useful to introduce the concepts such as the spin current and the 

spin torque to describe the coupling between the conduction electrons and the magnetic moments of ferromagnetic 

k> ' materials. Those are first proposed by Slonczewski [IjI based on a quantum-mechanical model. Theii, this interaction 

d i^ derived from the s-d coupling and manifests itself in magnetic multilayers 3J or bulk ferromagnets[l2- It also gives 

Ci ' a contribution as a current-induced force on a domain wall or an interaction between spin waves and itinerant 

electrons Q- The concepts are also used to deal with a variety of the structures such as the ferromagnet-normal 

metal- ferromagnet (FM-NM-FM) junctions [ij,[ia|j the ferromagnet-superconductor-ferromagnet junctions [igj and 

a trilayer FM-NM-FM contacting a normal metal lead or a superconductor lead '17', Ts^, and so on. Ap art from the 

quantum- mechanical method, the scattering matrix method is improved to cover this issue|l5l llTl Il8j| . Heide et al 

derived a set of coupled Landau-Lifshitz equations for the ferromagnetic layers by considering the nonequilibrium 

exchange interaction between layers'lOJ. Zhang et al dealt with the dynamic spin transfer torque and thermally 

assisted magnetization reversal by means of the Landau-Lifshitz-Gilbert equation 20, 21] in which a term describing 

the spin transfer torque is introduced. The nonequilibrium Green function is also applied to this question to consider 



the spin-flip scattering effect on the current-induced spin transfer torque 225- 

So far, the electrical current, the spin transfer effect and the spin current in magnetic multilayer systems are 
extensively studied under an dc bias voltage, and the investigation under an ac bias is still sparse. Whether an ac 
bias voltage applied to the magnetic tunnel junctions could induce some unusual properties (see e.g.J23), is still not 
clear. In this paper, we shall show the effect of an ac bias on the electrical current and the spin transfer torque in 
a magnetic tunnel junction. Since the switching of magnetic domains depends on the magnitude and the directions 
of the spin transfer torque, while the direction of the spin transfer torque is related to the direction of the electrical 
current, it is of interest to anticipate that the response of the system under an ac bias is different from that under 
an dc bias. We have found that the spin transfer torque is sharply enhanced by applying an ac electrical field with a 
particular amplitude and frequency. 

The rest of this paper is organized as follows. In Sec. II, a model is proposed, and the necessary formalism under 
an ac bias is established. In Sec. Ill, the tunneling current under an ac electrical field is derived and studied. In Sec. 
IV, the current-induced spin transfer torque under an ac electrical field is investigated. Finally, a brief summary is 
given in Sec. V. 

II. MODEL AND FORMALISM 

1. Model 

Let us consider a magnetic tunnel junction in which two ferromagnets (FM) which are stretched to infinite are 
separated by a thin insulator (I) . The molecular field in the left ferromagnet is assumed to align along z axis which is 
in the junction plane, while the orientation of the molecular field in the right ferromagnet deviates the z axis by an 
angle 6, which is along the z' axis (such that the frame x'oz' deviates xoz by an angle 9). The tunnel current fiows 
along the x axis which is perpendicular to the junction plane. The Hamiltonian reads 

H = Hl + Hr + Ht, (1) 

with 

Hl = ^eLWaLflfc-. (2) 

ka 

Hr = ^[(£qi?(t) - aNh cos 0)cl^Cqa ~ M2 sinOcl^Cq^], 

kqaa' 

where a^a and cj.^ are annihilation operators of conduction electrons with momentum k and spin a (= ±1) in 
the left and right ferromagnets, respectively. When the time-dependent bias voltage is applied, the single-particle 
energies in left and right ferromagnets, e^.^ (i), become time-dependent |24 . l25l | : e^^{t) = ekL(i) — crMi, ekL(^) =£fcL + 
AL{t)-eVo, Ai(i) = cVl cosc^q^, M, = ^^, s^R{t) =6^^ + Afl(t), Afl(t) = eVRCostoot. M^ = ^^, where Vo is 
the applied dc bias, Ai(i) and Aij(i) are from the applied ac bias, g is the Lande factor, fiB is the Bohr magneton, 
fT-L(R) is the magnitude of the molecular field of the left (right) ferromagnet, e^^, ^x is the single-particle dispersion 

of the left (right) FM electrode, T^f' denotes the spin and momentum dependent tunneling amplitude through the 
insulating barrier. Note that the spin-flip scattering is included in Ht when a' = a = —a. It is this term that violates 
the spin conservation in the tunneling process. 

2. Green Functions of Uncoupled leads 

First let us write down the lesser Green function and the retarded (advanced) Green function for the isolated leads 
which will be used subsequently. The single-particle energies of the isolated leads for spin up and down are splitting, 
i.e., Ef.'^ — EkL ^ (jMi{a ~ ±1, corresponding to t, i)- Following the standard procedure, it is not difficult to obtain 
the lesser Green function for the left FM lead 



—iujQ{mt—nt') 



(3) 



where /i(x) is the Fermi function of the left uncoupled lead, Jm(^^) is the m-th order of Bessel function. Throughout 
this paper, ?i = 1 is assumed. Similarly, we can write down the lesser Green function for the right lead. The retarded 
and advanced Green functions have the form of 

g;W(i,t') = Tj0(±tTi')VJ™( — )J„( — )e±^"°^"*-"*'^- (4) 

/ ^±^{el°-eVa){t-t') q \ 

where a — L,R. A further evaluation will make use of the double time Fourier transform as follows |2a 127] 

/OO /-OC 

dt, dt2F{ti,t2)exp[t{Eit,-E2t2)]. (5) 

-OO J — OO 

The above equations will be invoked in the following analysis. 

III. TUNNELING ELECTRICAL CURRENT UNDER AN AC BIAS 

In this section, we shall investigate the tunneling electrical current in FM-I-FM junctions under an ac electrical 
field. Following the method in Refs. |2J,|23j the tunneling current can be expressed as 

/(i)^^Re^T,v[flG<,(t,t)], (6) 

kq 

( T T \ 
where ft = TR, T = „ „ I with the elements Ti [i = 1, ...,4) of the tunneling matrix which are assumed to 

be independent of k and q for the sake of simplicity |2Sj, which is reasonable in the assumption of a wide-band limit 

(WBL)[2l|2l|23, and the elements Tz and Tg describe the effect of spin-flip scatterings [s^, R = ( *^°^ | ~ ^^V 

Tr^ stands for the trace of the matrix taking over the spin space, and G^ (t,i') is the lesser Green function in spin 
space defined as 



GiYit.t') Gi;'<(i,i') j' 



Gt,{t,t')^\ JU);' J!,< ', , (7) 



where G^°^ '^{t,t') = i{c^ (t')aka'it))- By using the noncquilibrium Green function technique |2J,[2^, we can get 
G^ {t,t') to the first order: 

where g]!^, g^ji, dkL-^ ^^"^ 9q^r ^^"^ retarded, advanced, lesser Green function of the left and the right ferromagnets, 
respectively. Then, substituting ((SJ into ((HJ, and making the double-time Fourier transform of Green functions, we 
have 



^W = T ^ j dEJ{t:')[fL{E") fn{E')]A{E',E'\9)cos{lcjot), 



(9) 



rnnm'ri' 



where J{Z;^') = J„(^)J„(£ii)J„,(£iii) j„,(£iii), E" = E - mc^o + eVo, E' = E - n'u^, I = m + m' - n - n' , 
K{E',E",e) - Tr„['[lT>R{E')n^TiL{E'% D^j^) = ( ^^^^^^ ^^0^^^ \ ^jti, DL(R)ui){e) = i?L(fl)(£ ± Afi(2)) the 

density of states (DOS) of the conduction electrons with spin up and down in the left (right) ferromagnet. In 
Eq. 0, the real part of the current is taken. The time average of the electrical current I{t) can be defined as 

{lit)) = — J \fy dtlL{t) with To the time interval between which the physical quantity is measured. Thus, the 

time-averaged tunneling current is given by 

javer ^ /jj-^^^ ^ f Y.mnm'n' ^im^n') = J2m7un'n' ^imn )IdciVo - muJo/e, n'uJo/e), I = 0; ^^q^ 

where liZ^^,) = J{Zn)Idc{Vo - m^o/e, n' u^ / e) , hc{eVo - muJo,n'LUo) = ^ J dE[hiE") - fn{E')]Tj{E', E"){1 

+y2{E )[n{E )COS0 + i'3[E jSmWJI, Fi - D^T(Ti'+T'|) + C,,i(T|+T|) ' ^2 - d„t+Chj' ^3 - D^,(Ti+ri)+D,,^(Ti+Ti)^ 

and r/(£;',£;") = [-DflT(^') + Dri{E')][Dl](E"){T^ + r|) + DLi{E"){Ti + T^)]. It can be seen that the average ac 
tunneling current in a magnetic tunnel junction is modulated via Bessel functions. This result is quite similar to Eq. 
(3.3) in Ref.|3l| in which the superconductor tunnel junctions are investigated. By analogy, the time dependence of 
the wave function for every electron state is modulated by J„( — £La.)g±™'^o*, Each single-particle level is modulated 
in terms of a probability J„( — i^^) and is displaced in energy by nhujQ. These displacements in energy contributing 
to the amplitude of the average tunneling current are equivalent to that dc voltages (Vb — muJo/e) applied across the 
junction with a probability J(JJ^„" ). In Ea. (|10|) . an explicit relation between the time averaged ac tunneling current 
and the dc current is given, where I = gives a nontrivial result. This result implies that only the elastic transmission 
through the tunnel barrier contributes to the average current, and the net number of photons absorbed from the ac 
field must be zero|23- So, every term in the summation of Eg. 1)10(1 describes the tunneling process that the electrons 
tunnel from the excited states eVo — muj absorbed m photons of the left ferromagnet to the excited states n'ujo absorbed 
n' photons of the right one. In this way, the summation is taken over all the excited states of the left and the right 
ferromagnets. The ac field gives a correction to the transition rate F/ by adding the product of Bessel functions, i.e. 
J{mn )i which describes the probability of electron population in the excited states. The tunneling current can be 
viewed as to happen between the photonic sidebands of the two ferromagnets. 

Next, let us present the numerical results of the time-dependent electric current. Before going on, we shall first 
give some presumptions. A parabolic dispersion for band electrons is assumed, on which is based that the DOS of 
conduction electrons are calculated. The Fermi energy and the molecular field will be taken as Ef = 1.295 eV and 
|hi| — |h2| — 0.90 eV, which are given in Ref.|33| for Fe. We note that the elements of the coupling matrix T2 and 
T3 mean the strength of spin- flip scatterings which were discussed in Refs. |23,|33- In order to focus our attention 
on the effect of an ac bias, we shall not consider the effect of spin-flip scatterings here for brevity and simplicity. 
So we take T2 = T3 — 0, and Ti — T4 = O.OleF. In addition, Vr = is assumed, then Vl = Vac- Under this 
assumption, I™'^'' = X)m •^m(^^)^rfc(^ — 'mujQ/e,6), where Jm{^^^) is the m-th Bessel function. When the ac 



LJQ 



bias is absent, we get GdciVo = 0,T = 0,6*) = dIdc/dVo = Go[l+ P2 \/P[+Pi cos{e - 0f)]^ at zero dc bias, zero 
temperature (T = 0) and 9, where Pi,P2, and P3 are similar to the formulas as given before, tan6'/ = P3/P1 and 
Gq = ^ [{T? + Ti)DL^ + {Ti + TI)Dli]{Dr^ + Dri). As a consequence, hc{T = 0,0) = Gdc(0, 0, 9) ■ Vq can be used 
as a scale to measure the ac current. 

The time evolution of the tunnel electrical current I{t) in response to an ac field is depicted in Fig. 1, where 
Vb = — O.iy and wq = 0.003 eV . As can be seen, when the dc bias is positive, the current can be negative. When 
the ac signal is small (e.g. Vac = O.OOIV), the current varies with time in a cosine manner, which appears to be 
proportional to the ac bias. In this case, the current response is similar to the dc case, i.e. larger the bias, larger 
the current. However, when the ac signal is stronger, things become different. There appear some resonant peaks, 
which can be regarded as to be resulted from the photon-assisted tunneling, and the tunnel current is still periodic 
with time. The inset of Fig.l shows the case at Vac — O.OIV^ but with a lower frequency for a comparison. It is found 
that the peaks split into several peaks in one period of time, suggesting that the external frequency can change the 
oscillating frequency of the tunnel current. 

The frequency dependence of the averaged tunneling current is shown in Fig. 2. Small oscillations of the current 
with frequency can be seen, which are caused by the summation of the m-th current /(^r^/) with different m, n, 
m', and n' owing to I{^n') oscillating with the frequency and having many peaks. The peaks of the m-th current 
correspond to the troughs of the (to -|- l)-th current. Furthermore, when the ac signal is small (e.g. Vac — O.OIF), 
the current first increases with small oscillations, and then, the current almost does not vary with the frequency like 
a pure resistance because the pure resistance should not vary with the frequency in the common sense. However, it 



can be observed that a large ac signal leads to a small averaged current also as those expressed in Fig. 3. The reason 
for this feature is that an ac bias is imposed on an dc bias, while the summation current decreases with increasing 
ac bias. A large Vac means a large argument of Bessel function, thus leading to a strong modulation to the current, 
which makes the current approaching to zero. If Vac is fixed, a larger loq gives a larger averaged current, suggesting 
that the current flows easily in this system under a higher frequency. This character is consistent with the classical 
feature of systems that metallic leads are separated by an insulator. Here, we can consider two limiting cases. When 
ujo -^ 0, the argument of the Bessel function xb = eVac/uJo — > oo, then J^{xb) ~ :;;J— cos'^{xb — f m — j) ^ 0, which 
suggests that the current cannot pass easily through the system. When Vac —^ and the frequency becomes larger, 
the argument of the Bessel function xb ^ 0, then Jm^oi^B ^ 0) '^ and Jo(0) ^ 1, which suggests that the m = 
term is dominant. This case shows that the current can pass easily through the barrier under higher frequencies. It 
is the character of a capacitance in the usual sense. 

Let us define TMR — {Jp — Jap)/ Jap by using the real part of the current. We have investigated the response of 
TMR to an ac bias. It is found that TMR almost does not alter with Vac and frequency wqj with the varying range 
about in 0.01% ~ 0.1%. To understand this result, we would like to remark that TMR is mainly contributed by the 
spin-dependent scatterings. When electrons from one ferromagnet enter into another whose magnetization deviates 
an angle to the first one, spin up and down electrons bear different potentials. While the ac field provides the same 
modulation of quasiparticle levels of the spin up and spin down electrons shifted by nfwjQ^ leading to a change of the 
magnitude of the current, it affects less spin-dependent scatterings. If we impose a time-dependent magnetic field 
on both electrodes, it is conceivable that TMR would be remarkably influenced by such an ac magnetic field. Work 
towards this direction is now in progress. 

IV. CURRENT-INDUCED SPIN TRANSFER TORQUE UNDER AN AC BIAS 

As is known, the relative orientation of magnetizations on both electrodes can affect considerably the magnitude 
of the electrical current flowing through the magnetic tunnel junctions, and meanwhile, the spin-polarized electrical 
current can also switch the direction of magnetization of electrodes. This latter effect comes from an indirect inter- 
action between ferromagnets, which is caused by the so-coined current-induced spin-transfer torque (CISTT)[3|, |j]. 
The physical meaning of the CISTT can be understood as follows. When the electrons in the first ferromagnet tunnel 
through the barrier and enter into the right ferromagnet, the incoming polarized electrons will precess, eventually to 
align with the magnetization direction at an angle 6 in the second ferromagnet. In this process, there should be a 
difference of spin angular momenta between the incoming and outgoing spins in the second ferromagnet. The missing 
spin angular momentum must be absorbed by the local moments, thereby generating a torque which exerts on the 
moments of the second ferromagnetic layer. This kind of torque reflects actually the spin angular momentum transfer 
from the first ferromagnet to the second by the polarized current. In Ref. [23, we have discussed the CISTT by 
means of the nonequilibrium Green function technique under an dc bias. By generalizing our treatment, we now 
consider the case under a simultaneous application of both ac and dc biases. The CISTT is related not only to the 
available energy levels of both ferromagnets, but also to the direction of the current. In the ac case, the available 
energy levels of ferromagnets will be shifted by absorbing or emitting photons. The tunneling can be viewed as to 
occur between the photonic sidebands of the left and the right ferromagnets, and the magnitude of the current can 
be tuned. The ac field can change the current directions by frequency 2a;o- It may be expected that the ac field may 
impose considerable effect on the CISTT. The total spin of the right ferromagnet can be expressed as|2^ 

«W = I E 4c..(R"V,)t a (R-V.), (11) 

where o is the Pauli matrices and X^i{v) is spin states, which is written down in the xyz coordinate frame while the 
spins S2 are quantized in the x'y' z' frame, and the rotation matrix R is the same as before. Since 



*L,i?/- 



Ic'SL,R X (SL X Sfl), (12) 



where s'l,r are unit vectors Si = Si/si and hsi represents the respective total spin momenta per unit area of the 
ferromagnets Q, we know that the direction of the spin transfer torque Sp is just along the x' direction in the x'y' z' 
coordinate frame. From Eq. (fTT|) we obtain S2(t) = f X] {c\^Cka cos 9 — ac\^Cka sva 9) = S2a;'o cos0 — S2z'o sin0, where 

k(T 

S2x'o and S2z'o are x'- and z'-components of the total spins in the x'y' z' coordinate frame in which the spins S2 are 



quantized. So the CISTT can be obtained 22^ 

r^-(i) = -cos^Re ^ Tr,[G<^(t,t) a^ T^] +sin0Re ^ Tr,[G<g{t,t) a^ T^], (13) 

kq kq 

where cri= in)' '^3^ n i ) ^^^ Pauh matrices, and G"^ {t,t') is the lesser Green function in spin space 

defined as above. By using the noncquihbrium Green function technique [2jj (a similar framework employed in 
Ref . |23| ) . we can get the torque to the first order of the Green function 

G,< (t,i') = y'dii[g^i(i,ii)TRg<,(ti,i')RUg,<L(t,ii)TRg^«(ii,i')R^], (14) 

where gJJ^, Sgji, SkL' ^^^ S^fl ^^''^ retarded, advanced, lesser Green function of the isolated left and the right ferro- 
magnets, respectively. After some algebra, we obtain the Fourier transform of the CISTT 



T 



R^i,,\ - o^2 Y^ 



mnm'n' 



M = 2^^ > dEJ{2-)[f^{E")~fn{E')]nE',E" 



■[Pi{E")sme~ PaiE") cos 0]S{uj + Iujq), 

where T{E', E") = [Dii]{E') + DRi{E')][DL]iE"){Tl+T^)+DLi{E")(Ti +Ti)]. Here we may note that the direction 
of the spin torque is related not only to whether the applied dc bias is positive or negativeja. Isl rTlla. Il(]| (namely, the 
direction of the electrical current in a steady state) , but also to the frequency of the external ac bias and the Bessel 

la 
function's order. The time average of CISTT may be defined as (t^^) — ^ J %^ dtT^^{t). Then, we get 



R^ ^ (^B.^\ ^ [ ^j:„.nm'n'IdEJ{±-')[fL{E") - f r{E')]V{E' , E")P{E" , 6) , I = 0; , . 



T 
aver 



where I = m + m' - 7i - n' , E' , E", J{Zn')^ P{E", 0) = [Pi{E") smO - PsiE") cose*], and r(^', E") are defined as 
above. 

The 6 dependence of the spin-transfer torque is similar to that in Refs. |a. Ha . UtI lia | . In a colinear case (0 = or tt) 
and without spin-flip scatterings, the CISTT disappears, even though an ac bias is present. At 9 = ^, the CISTT tends 
to a maximum in the presence of no-flip of spins 22J . In this case, loq, Vac, and Vq can affect the maximum of the torques. 
In the following discussions, we suppose that Ti = T4, T2 = Ts = 0, Vr = and Vl — Vac- We shall use Tdc{T = 0, 6*) = 
GS,(O,O,0) • K) as a scale, where G5,(O,O,0) = eTTTf[DR^{Ef) + DRi{Ef)][DL^{Ef) + DLi{Ef)]Pi{Ef) sm0 with 
P^{Ef) = [DL^iEf)~DLi{Ej)]/[DL^iEf)+DLiiEf)]. 

First, in order to observe the effect of an ac bias on the CISTT, we present Vac dependence of the time-dependent 
spin-transfer torque at different frequency wq in Fig. 4. It is seen that the amplitude and oscillating frequency of the 
torque are modulated. As the direction of the spin-transfer torque is related to the direction of the electrical current, 
it is changed continuously with the ac current. Note that in Fig. 4, an dc voltage is applied simultaneously together 
with the ac bias. When Vac is small, the oscillating amplitude becomes large, and when Vac becomes larger, the torque 
is strongly suppressed. Besides, one may see that a larger loq gives a small oscillating frequency of the torque with 
Vac- The frequency (wq) dependence of the time-dependent CISTT at different amplitudes of the ac bias is shown in 
Fig. 5. It is seen that the torque oscillates with the external frequency, and it appears that the signal of the CISTT is 
strong in some regime of wqi and is almost vanishing in other regimes. This feature becomes more evident at a larger 

Vac- 

The time-averaged spin torque (t^^^\ is more interesting. In Fig. 6, the Vac dependence of the time-averaged 
CISTT is plotted under different frequencies of the ac bias. One may find that the time-averaged CISTT increases 
slowly with small Vac, and then decreases rapidly with small oscillations. The CISTT approaches to zero when Vac 
becomes larger. As shown above, it seems that a large amplitude of the ac bias field may suppress the spin-transfer 
torque. When Vac is decreasing, (t^^) /''dc(O) approaches asymptotically to a single curve for different frequency ljJq. 
At a given wq, (r^^) /T'dc(O) has a maximum at a specific Vac- In this case, the time-averaged CISTT, <^t^^) (note 
that we use t^^^^^ to denote it in figures herefater), in the simultaneous presence of an ac bias and a dc bias, can be 
expressed by that under only a dc bias: 

(r«-)= Y. JCl"')r,?(^o-nWe,mc.o/e), (16) 

mnm'n' 



where the subscript dc represents the quantity under only a do bias. It is clear that the Bessel function modulates 
the amplitude of the torque. From Eq. ((T^ . one may see that a large dc current generates a large torque. In the 
case under an ac bias, the electrical current might have the similar character. Since the tunneling is viewed as to take 
place between the modulated levels of ferromagnets, the modulations lead to different transmissions, thus enabling 
the spin-transfer torque to exhibit various features. This character also manifests itself in Fig. 7, in which the time- 
averaged spin-transfer torque as a function of the external frequency luq for different Vac is shown. It can be found 
that small frequencies affect the time-averaged CISTT dramatically, manifested explicitly by an approximately linear 
relation between the CISTT and small loq. In comparison to Fig. 2, it may be concluded that (T^^(t)) is proportional 
to the time-averaged electrical current, say {Ie{t)), and the CISTT can be still induced by the ac electrical current. 
Note that Eq. ((T^ was proposed under an dc biasQ. Under a small ac signal, the torque almost does not vary in 
higher frequencies, which is saturated with the magnitude larger than that under only a dc bias. 

In above, we have applied simultaneously both the dc and ac biases to the tunnel junction. In the presence of only 
an ac bias, the time-averaged CISTT as a function of Vac is given in Fig. 8. It is seen that the time-averaged spin 
transfer torque first increases sharply and then decreases rapidly with increasing Vac- In other words, the CISTT 
takes its larger values in a narrow regime of Vac- The higher the external frequency ujq, the larger the maximum of the 
torque. This is because the spin transfer torque is related to the subtraction of two Fermi functions, and the latter is 
very small when only an ac bias is present. The principal contribution to the result comes from a narrow regime of 
Vac- The inset of Fig. 8 gives a plot of V^^ versus luq, where VJ^ is the specific value of the amplitude of the ac bias 
at which the peaks of the torque appear. It is to note that VJ'c is almost proportional to the frequency luq. Compared 
to Fig. 6, where the dc and ac biases are present simultaneously, the peaks in Fig. 8 are more sharp. These peaks 
can be regarded as the evidence of the photon-assisted enhancements of the CISTT. In Fig. 6, the averaged spin 
transfer torque increases slightly with increasing Vac, and tends to a maximum, then decreasing rapidly. If the dc 
bias disappear slowly, the spin transfer torque tends to zero when Vac tends to zero, while the situation in Fig. 6 is 
not. This reveals that the CISTT exhibits different behaviors under different biases. Our results show that an ac bias 
could enhance the CISTT at particular values of the amplitude and frequency. 

V. DISCUSSION AND SUMMARY 

We have presented the calculation of the response of the electrical tunnel current and the current-induced spin 
transfer torque to ac and dc biases in FM-I-FM tunnel junctions by means of the time-dependent nonequilibrium 
Green function technique. In general, when a spin-polarized current is injected into a ferromagnet, the polarized 
electrons may experience the fields such as the external field, the anisotropy field, the dipolar-dipolar interaction, and 
the demagnetization field[i^|33, leading to the precession of the polarized electrons around these fields. To conserve 
the spin angular momentum, the magnetization also precesses around the polarized direction of the injected electrons, 
giving rise to a spin transfer which is quite different from a magnetic field induced by the current. The ways to realize 
the switching of magnetization can be either by a magnetic field, or by a magnetic field induced by the current, or 
by the spin transfer effect ,35i]. The way through the mechanism of spin transfer offers a convenient and fast choice to 
switch magnetization, which could be observed from other works and our present study. 

In summary, we have investigated the spin-dependent tunneling in the presence of an ac bias applied to a FM-I-FM 
system by considering the ac tunneling current and the ac CISTT. We have formulated the time-averaged current 
(time-averaged CISTT) in the form of a summation of dc current (dc CISTT) multiplied by products of Bessel 
functions with the energy levels shifted by mTiujQ. The tunneling current can be viewed as to happen between the 
photonic sidebands of the two ferromagnets. Our calculation shows that low- frequency ac field suppresses the current, 
and the electrons may more easily tunnel through the barrier under a high-frequency ac field as the response for a 
capacitance in a classical case. It is found that the TMR almost does not vary with an ac bias, which suggests that 
the ac electric field contributes less to the spin-dependent scatterings. The current-induced spin transfer torque under 
an ac bias has also been investigated. It has been shown that an ac bias may overall suppress the spin transfer torque, 
but in a narrow regime of the ac bias, the CISTT is greatly enhanced, characterized by a sharp peak which can be 
viewed as the photon-assisted enhancement. It has been found that the particular amplitude of the ac bias at which 
the CISTT shows a peak has a linear relation with the frequency of the ac bias approximately. As a consequence, 
people could adjust the proper values of the amplitude and the frequency of the ac bias in accordance with such a 
linear relation to realize the enhancement of the CISTT, which gives another possible option to enhance the CISTT 
externally. In addition to using an dc current, our result might give an alternative hint to control the local moments 
of ferromagnets, and also to control the resistance states. 
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FIGURE CAPTIONS 

Fig. 1 Time dependence of the tunneling current, where the effective masses of the left and the right ferromagnets 
are taken as unity, the molecular fields are assumed to be O.QeF, the Fermi energy is taken as 1.295eT^ which are taken 
from Ref.|33| for Fe, and Ti = T4 = 0.01 eV, T2 = T3 = 0, the other parameters are assumed as ?i = 1, Vq = —0.1 V, 
ujo = 0.003 eV, 9 = tt/S and the temperature is at lOOK. 

Fig. 2 The frequency dependence of the time-averaged tunneling current. The remaining parameters are assumed 
to be the same as those in Fig. 1. 

Fig. 3 The time-averaged tunneling current as a function of the amplitude of the ac bias. Vac- The remaining 
parameters are assumed to be the same as those in Fig. 1. 

Fig. 4 Time-dependent spin-transfer torque as a function of Vac under different frequencies, where t — 15 (lO^'^ns), 
Vq = —0.05V, and the other parameters are assumed to be the same as those in Fig. 1. 

Fig. 5 Time-dependent spin-transfer torque as a function of frequency uq under different Vac, where Vq = — O.IV, 
and the other parameters are the same as those in Fig. 4. 

Fig. 6 Vac-dependence of the time-averaged spin transfer torque under different frequencies luq, where the other 
parameters are the same as those in Fig. 1. 

Fig. 7 Time-averaged spin transfer torque as a function of luq for different Vac, where the other parameters are the 
same as those in Fig. 1. 

Fig. 8 Time-averaged spin transfer torque versus Vac for different frequencies when only the ac bias is present 
(Vq = 0). The inset is the particular amphtude, V^^, at which the CISTT is peaked, versus frequency ujq. The 
parameters are the same as those in Fig. 1. 
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FIG. 1: Time dependence of the tunneling current, where the effective masses of the left and the right ferromagnets are taken 
as unity, the molecular fields are assumed to be 0.9el/, the Fermi energy is taken as 1.295eV which are taken from Ref.|33l for 
Fe, and Ti = T4 = 0.01 eV, T2 = Ts = 0, the other parameters are assumed as R = 1, Vb = -0.1 V, cjq = 0.003 eV, 9 = tt/S 
and the temperature is at WOK. 
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FIG. 2: The frequency dependence of the time-averaged tunneling current. The remaining parameters are assumed to be the 
same as those in Fig.l. 



11 



S- 0,6 




FIG. 3: The time-averaged tunneling current as a function of the amplitude of the ac bias, Vac- The remaining parameters are 
assumed to be the same as those in Fig. 1. 
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FIG. 4: Time-dependent spin-transfer torque as a function of Vac under different frequencies, where t = 15 (10 '''ns), Vo 
—0.05V, and the other parameters are assumed to be the same as those in Fig. 1. 
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FIG. 5: Time-dependent spin-transfer torque as a function of frequency loq under different Vac, where Vq 
other parameters are the same as those in Fig. 4. 
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FIG. 6: Vac-dependence of the time-averaged spin transfer torque under different frequencies too, where the other parameters 
are the same as those in Fig. 1. 
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FIG. 7: Time-averaged spin transfer torque as a function of uo for different Vac, where tfie otfier parameters are tfie same as 
tliose in Fig. 1. 
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FIG. 8: Time-averaged spin transfer torque versus Vac for different frequencies when only the ac bias is present (Vo = 0). The 
inset is the particular amplitude, V^c, at which the CISTT is peaked, versus frequency ojo- The parameters are the same as 
those in Fig. 1. 



